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Hyperspherical Partial Wave Theory has been applied to calculate T-matrix elements and Single 
Diflterential Cross-Section (SDCS) results for electron-hydrogen ionization process within Temkin- 
Poet model potential. We considered three different values of step length to compute the radial part 
of final state wave function. Numerical outcomes show that T-matrix elements and SDCS values 
depend on the step length h. Here, we have presented T-matrix elements and the corresponding 
SDCS resuhs for 0.0075 a.u., 0.009 a.u. and 0.01 a.u. values of h and for 27.2eV, 40.8ev and 54.4eV 
impact energies. With the help of the calculated data for three different step lengths, we have been 
able to evaluate a two-term error function depending on the step length h. Finally, two-term error 
corrected T-matrix elements and the corresponding SDCS values have been computed. We fitted 
our two-term error corrected SDCS results by a suitable curve and compared with the benchmark 
results of Jones et al. [Phys. Rev. A, 66, 032717 (2002)]. Our fitted curves agree very well with 
the calculated results of Jones et al. and two-term error corrected SDCS results somewhere agree 
with the benchmark results. Two-term error corrected SDCS results are significantly better than 
the calculated SDCS results of different step lengths. 

PACS numbers: 34.80.Dp, 34.10.+X, 34.50.Fa 
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I. INTRODUCTION 



The electron-impact ionization of hydrogen probes the 
correlated quantal dynamics of two electrons moving in 
the long-range Coulomb field of a third body. As such 
it remains one of the most fundamental and interesting 
problems in nonrelativistic quantum mechanics. There 
are many attempts for a complete solution but all of 
these face enormous difficulties and have only limited 
success. Among these the most successful attempts are 
the method of Convergent Close-Coupling (CCC) and 
Exterior Complex Scaling (ECS) . Another promising ap- 
proach for the electron-hydrogen atom ionization prob- 
lem is the Hyperspherical Partial Wave (HPW) approach. 
After the successful applications of HPW theory to com- 
pute triple differential equal-energy-sharing cross-section 
results P, 0, H, I3j S @| J we aspire to calculate SDCS re- 
sults. Before considering the full electron- hydrogen ion- 
ization problem, here, we consider Coulomb three-body 
system within Temkin-Poet (TP) model 0,[i]. The TP 
model of electron-hydrogen collision is now widely con- 
sidered as an ideal testing ground for the improvement 
of general methods intended for full Coulomb three-body 
problem. In this context, the calculated SDCS results of 
other theories for TP model potential are praiseworthy. 
Among these the attempt of Jones et al. ^, Jii] is remark- 
able, they obtained benchmark results. They have devel- 
oped a variable-spacing finite-difference algorithm that 
rapidly propagates the general solution of Schrodinger 
equation to large distances, originally used by Poet [llj 
to solve TP model. The ECS calculation is generally 
in good agreement with the benchmark results of Jones 



et al. except at the extreme asymmetric energy shar- 
ing [13. The calculated singlet SDCS curves of CCC 
method are wavy, Bray considered a smooth curve by 
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educated guess [12|. The CCC results agree nicely with 
benchmark results of Jones et al. only for the triplet 
case (generally, CCC does not yield convergent ampli- 
tude for the triplet case, except for total angular mo- 
mentum zero). We also note the work of Miyashita et 
al. They have presented SDCS for total energy of 

4Ry, 2Ry and O.lRy using two different methods. One 
produces an asymmetric energy distribution similar to 
that of CCC while the other gives a symmetric distri- 
bution. Both contain oscillations. It should be noted 
that recently, we have used HPW approach to calculate 
SDCS results for full electron-hydrogen-ionization prob- 
lem at 60eV incident energy (15.] . The resultant curve 
was wavy and calculated cross-section results are irrele- 
vant at extreme energy sharing. We had fitted our cal- 
culated SDCS data by a fourth order parabola and com- 
pared with the experimental values of Shyn [l^ . Our fit- 
ted curve agrees excellently with experimental results. In 
this article we present the SDCS results for TP model us- 
ing HPW method with two-term error correction. Here, 
we introduce a procedure to calculate error function. The 
results are obtained for intermediate (27.2eV, 40.8eV and 
54.4eV) energies. We have calculated T-matrix elements 
and the corresponding SDCS data for three different val- 
ues of step length h (0.0075 a.u., 0.009 a.u. and 0.01 a.u.), 
use to calculate radial part of final state wave function 
numerically. Numerical observation shows that the T- 
matrix elements depend on h. Using the data for various 
step lengths, we calculated two-term error function, de- 
pends on h. Finally, two-term error corrected SDCS val- 
ues were computed. The nature of error corrected SDCS 
undulating curves suggests a fit, with a proper function. 
HPW method for TP model is reproduced in Sec. II, pro- 
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cedure of calculation is presented in Sec. Ill, two-term 
error correction process is given in Sec. IV, results are 
presented in Sec. V with a short discussion, and some 
concluding remarks arc found in Sec. VI. Atomic units 
are used throughout this paper except where otherwise 
noted. 



where 
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FIG. 1: (Color online) The values of T^„i (zero indicates sin- 
glet) for three different step lengths at 27.2eV incident elec- 
tron energy. Square points for h — 0.0075, diamond points for 
h = 0.009 and pentagon points for h = 0.01. Hexagon points 
represent the values of y^'J^' at the same energy. 



II. THEORY 

The T-matrix element, we use in cross-section calcula- 
tion, is given by 



(1) 



In this expression $i is the unperturbed initial channel 
wave function, satisfying certain exact boundary condi- 
tion at large distance and Vi is the corresponding pertur- 
bation potential. Here, '^^fj is the symmetrized scatter- 
ing state (see Newton for definition) . For information 
regarding electron-hydrogen-ionization within TP model 
potential, one may solve the corresponding Schrodinger 
equation. We start by writing the Schrodinger equation 
for the full electron-hydrogen ionization problem 



1 



Vi2 = ^ ^. (3) 

\ri - r2\ 

To calculate the final channel symmetrized continuum 
state '^'f'J we use hyperspherical coordinate R = 
V'^i+'^h = arctan{r2/ri), r'l = (6*1, 0i), r2 = (6*2, 
and uj = (a,fi,r2). Also we set P = a/pi+pI' 
ao = arctan{p2/pi), pi = {0p^,(f)p^), P2 = (^p2)<^p2) and 
^0 — (ao,Pi,P2) where fl and pi {i — 1, 2) are the coor- 
dinates and momenta of i th charged particles, '^^fj is 
then expanded in symmetrized hyperspherical harmon- 
ics [H that are functions of five angular variables and 
li,l2,n,L,M, which are, respectively, the angular mo- 
menta of two electrons, the order of the Jacobi polyno- 
mial in hyperspherical harmonics, the total angular mo- 
mentum and its projection. For a given symmetry s (s 
= for singlet and s = 1 for triplet), we decompose the 
final state as 



(-) 



E 



.5/2 



(4) 



where /i is the composite index {li,l2,n, L, M) and p = 
PR and 4>1i_{io) are orthogonal functions that are product 
of Jacobi polynomial P;";^ and coupled angular momen- 
tum eigenfunction Y^^/^^ {fi,r2) F^{p) then satisfy 
the infinite coupled differential equations 



1- 



F^,{p) = 0. (5) 
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E^'jJ (2) FIG. 2: (Color online) Same as in Fig. 1 but for triplet case. 
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In the above expression 



FIG. 3: (Color online) Same as in Fig. 1 but for 40.8eV 
incident electron energy. 



Here af,,,/ are the matrix elements of the full three-body 
interaction potential and vx = A + 3/2 where X = 2n + 
h+h- 

For the cusp model (or TP model) the V12 term, de- 
rived from the first term of the partial-wave expansion of 
the electron-electron potential, is given by 



V12 



1 



1 



r> max{ri,r2) 



(6) 



The TP model calculated in this article is simplification 
of our earlier calculated full electron-hydrogen problem, 
and we only consider the case where all angular momenta 
are zero. Retaining only zero angular momentum terms 
we have 



(-) 



(7) 



where 



The expression of hyper- 



spherical harmonics where all angular momenta are zero 
is given by 

0(L=Z,=Z.=O),„M = + (-l)^+"}Po"o(«)^00°(''"l,''2). 

(8) 

The radial functions (p) satisfy an infinite coupled set 
of equations 



dp' 



F,1{p) + Y,^F:Ap)=0. (9) 



and 



C 



cosa sina max{cosa, sina) 



(10) 



(11) 



Finally, one obtains the T-matrix element in the form 
(for details see Eqn. (25) of Ref. [l|) 



(12) 



The modulus square of the T-matrix element, which is 
used to calculate differential cross-section, is then given 
by 



III. PRESENT CALCULATION 

In our present calculation, n, the degree of Jacobi poly- 
nomial, was varied from to 11. We considered n = 0, 2, 
4, 6, 8, 10 for calculating singlet SDCS results and n = 
1, 3, 5, 7, 9, 11 for computing triplet SDCS values p^ . 
The main numerical task is to calculate the radial func- 
tions F^p) over a wide domain [0,oo). As earlier we 
divide the whole solution interval [0, 00) into three subin- 
tervals [0, A], (A, Rq] and [i?o, 00), where A has the value 
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FIG. 4: (Color online) Same as in Fig. 3 but for triplet case. 
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FIG. 5: (Color online) Same as in Fig. 1 but for 54.4eV 
incident electron energy. 



of a few atomic unit and Rq is the asymptotic matching 
parameter. Rq is needed to be such that Rq ^ 1/VE, 
where E is the energy in the final channel 1]. Thus for 
energies of 27.2eV, 40.8eV and 54.4eV this range param- 
eter i?o may be chosen greater than the values 5000 a.u., 
3000 a.u. and 2500 a.u., respectively. We have chosen Rq 
around these values in our calculations. For [i?o,oo) we 
have simply analytic solution We applied a seven- 
point finite difference scheme [3| for solution in the in- 
terval [0, A] with step length h. Now for the difference 
equations we divided the domain [0, A] into 100 subinter- 
vals of length h and A = lOO/i. Solution over (A, Rq) is 
very simple. Because of the simple structure of equation 
(9) a Taylor series expansion method with step length 
2h works nicely. Presently, we considered three different 
values of step length h, these are 0.0075 a.u., 0.009 a.u. 
and 0.01 a.u., respectively. Finally, we calculated T^„, 
and SDCS results for three different step lengths. 



IV. TWO-TERM ERROR CORRECTION 

In the previous section, we reproduced the values of 
^rm' ^'^^ three different step lengths and observed that 
-^nn' varied with h. Now, we can consider a relation 
between T^„, (h) with the error term E^^^, (h) as 

T^n'ih)^TZ>+E:,^,{h) (14) 

where T^*, are the converged results with respect to the 
step length h. Since in our seven-point finite difference 
scheme the error term is 

Kh^f ^XO H where K IS a 



constant and ^ is a linear function of h. The error term of 
T^^ni{h) calculation is Ch^f'-^\(,) where C is a constant. 
Instead of Kh^f'-^^£_), we can write the error term of 
seven-point finite difference scheme as 

K,h'f^'\R,r.)+K,h'y^''^\Rr.)+K,h''f(''\0, 

for a certain grid point Using the above expression, 
we can easily formulate, 

Kn'W = ^L'/i' + B-L,h^" + G-i,h^'f^''Ho (15) 



where _B^„, and G'^„, are independent of h. Consid- 
ering first two terms, we get the expression of two-term 
error function for T^„i{h) elements 
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(16) 



Corresponding two-term error corrected T„„, (h) ele- 
ments satisfy the equation 



(17) 



In the present context, we have considered step lengths 
of three different values /ii, /12 and /13. Therefore, from 
the equation (17) we have, 

K^nhh.) - E2}{h,) = T^„,{h.) - T^„,{h,) (18) 

for i, j = 1, 2, 3 and i^j. The coefficients of and 
in the expression (16) are given by 
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FIG. 6: (Color online) Same as in Fig. 5 but for triplet case. 
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FIG. 7: (Color online) Singlet SDCS {nal/Ry) vs the en- 
ergy fraction Et/E for three different step lengths and for 
r°*S'^' elements at 27.2eV impact electron energy. Continu- 
ous curve, elements; dashed curve, for h = 0.0075 a.u.; 
dash-dotted curve, for h = 0.009 a.u.; dash-double dotted 
curve, for h = 0.01 a.u.. 



{hf - hf){T^,Ah2) - TL.{K)} ^ {hf - hf){T:^M - T^h,)} 
{hl^hl){hf-hf)-{hl-hl){hf-hf) 
{hi - hl){T,\^,{h^) - T^^Ahi)} - {hi - h\){T,l^,{h^) ~ r,t„,(/.2)} 
{hl-h\){hf-hf)-{hl-hl){hf-hY>) 



Here we have considered ft,i=0.0075, ft-2=0.009 and duce to 
/i3=0.01 so the above expression of v4f„j, and re- 



= 0.05229064077{T^„,(/i2) - n^.{hi)} - 0.023472188{r4,(/i3) - n^.{h2)} 
B^, = -0.01143107936{r4,(/i2) - n^>{hi)} + 0.006630530581{r4,(/i3) - n,,{h2)}. (20) 



After calculating the T,^*; ' elements, we have calculated 
the corresponding two-term error corrected SDCS re- 
sults. 



V. RESULTS AND DISCUSSION 



As we discussed in the section III, we have considered 
six different values of the degree of Jacobi polynomial. 
There are total 36 pairs of (n, n') in the calculation of 
T^„, . We have labeled those pairs by an integer variable 
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FIG. 8: (Color online) Triplet SDCS {-Kal/Ry) vs the en- 
ergy fraction Et/E for three different step lengths and for 



FIG. 10: (Color online) Same as Fig. 8 for 40.8eV. 
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FIG. 9: (Color online) Same as Fig. 7 for 40.8eV. 
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FIG. 11: (Color online) Same as Fig. 7 for 54.4eV. 



P, varied from 1 to 36. The values of T^^' (zero indi- 
cates singlet) for three different step lengths and T^*^^ 
are presented in Fig. 1 for 27.2eV energy, in Fig. 2 for 
40.8eV energy and in Fig. 3 for 54.4eV energy. In the 
Figs. 4, 5, and 6 we have presented the values of T^„, 
(one indicates triplet) for three different step lengths and 



7 




0.2 0.4 0.6 0.8 1 

energy fraction 



FIG. 12: (Color online) Same as Fig. 8 for 54.4eV. 
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FIG. 13: (Color online) Singlet SDCS [iral/Ry) vs the energy 
fraction Ei,/E for incident energy 27.2eV. Continuous curve, 
fitted function; dashed curve, present results corresponding 
rp0*(2) gjgjjjgjj^g. dotted curve, calculated results of Jones et 

ai. [ig]. 



r„^*P^ for energies of 27.2eV, 40.8eV and 54.4eV respec- 
tively. Figures show that the magnitudes of T^„, are 
diminished with the decreasing of the step length. As 

shown in the figures, the magnitudes of T^„, are low- 
est than the values of T^„, for various step length. The 
curves were drowned joining the points square for h = 
0.0075, diamond for h = 0.009, pentagon for h = 0.01 
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FIG. 14: (Color onhne) Triplet SDCS {ival/Ry) vs the energy 
fraction Eb/E for incident energy 27.2eV. Continuous curve, 
fitted function; dashed curve, present results corresponding 
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FIG. 15: (Color online) Same as Fig. 13 for 40.8eV. 



and hexagon for T^J) , show that in the figures, com- 
paratively similar. In our previous calculation primarily 
to calculate Double Differential Cross-Section (DDCS), 
SDCS results and somewhere Triple Differential Cross- 
Section (TDCS) results, we had established good quali- 
tative results. There were significant discrepancies in the 
magnitude for extreme asymmetric energies. These types 
of phenomena were happened due to such kind of behav- 
ior of T^jj, elements, depend tremendously on h. At that 
time, we drew full electron-hydrogen problem, and it was 
difficult to envisage the convergence analysis with respect 
to step length. These figures also show that the calcu- 
lation of Tf^^, elements is stable concerning h. The two- 
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FIG. 16: (Color online) Same as Fig. 14 for 40.8eV. 
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FIG. 17; (Color online) Same as Fig. 13 for 54.4eV. 



term error corrected elements T^^, are almost less than 
T^„, (0.0075) and in few cases equal with r^„, (0.0075). 
This implies that tow-term error correction procedure 
will be fruitful. In the section, we shall show that the 
SDCS results for r„;r' elements are significantly better 
than the SDCS results for h = 0.0075 and other values. 
In Figs. 7-12, we have compared our calculated SDCS 
results for three different step lengths and for T^j*!^^ (s — 
for singlet and s = 1 for triplet) elements. As shown in 
figures, the SDCS curves are less corrugated and smaller 
magnitude with the decreasing of step length. The calcu- 
lated SDCS results for T;;r' elements are smallest in size 
and least undulating comparison with the SDCS values 
for three different step lengths. In the case of singlet for 
27.2eV and 54.4eV energies, the SDCS results for t"*^^ 
elements are significantly different to that of for three 
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FIG. 18: (Color online) Same as Fig. 14 for 54.4eV. 

different values of h. At 54.4eV energy, there is an irrele- 
vant peak at equal energy sharing case for h = 0.009 and 
0.01. The peak abolished for h = 0.0075 and reduced to 
a deep for T„„; elements. For 40.8eV energy, the shape 
of the curves is approximately same only difference in 
magnitude. Same things happened in the case of triplet, 
the scale of curves reduced and wavy nature abolished, 
as well as singlet. With the change of step lengths, the 
nature of curves change rapidly, its magnitude reduced 
and wavy nature abolished swiftly with the modification 
of the step length, in decreasing order. Magnitude of the 
curves reduced significantly at the extreme asymptotic 
region. 



VI. COMPARE WITH BENCHMARK RESULTS 

In this section, we present two-term error corrected re- 
sults and fitted curves corresponding these values along 
with the benchmark results of Jones. The oscillating na- 
ture of two-term error corrected curves suggests a fit, 
with a proper function, symmetry about E/2 (E is the en- 
ergy in the final channel). We looked at the Unear-linear 
function for singlet SDCS values {y — a + bx + c\x — d\ 
where x is the energy of the secondary electron and y 
is the corresponding singlet SDCS values) and a maxi- 
mum six degree polynomial for triplet SDCS data (y — 
a + bx + cx^ -f dx^ -I- ex'^ -I- Jx^ -V gx^ where x is the en- 
ergy of the secondary electron and y is the correspond- 
ing triplet SDCS values). For 27.2eV energy, a four de- 
gree polynomial proved sufficient for curve fitting. First, 
we neglected broader data (maximum eight data out of 
eighty) from the extreme asymptotic region, irrelevant 
with other data, fitted a function for the rest of the val- 
ues and drew the fitted curve for the entire energy do- 
main. In Tables 1 and 2, we have presented the coef- 
ficients of the fitted curves for three different energies 
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TABLE L Coefficients for fitted curve at various impact energies and tripiet case. 



Energy 


a 


b 


c 


d 


e 


f 


g 


27.2eV 


0.042012 


-0.650385 


3.8148108 


-10.0556766 


10.054884 






40.8eV 


0.046054 


-0.49753 


2.90446 


-9.505798 


16.482738 


-14.075808 


4.691936 


54.4eV 


0.0563445 


-0.420831 


1.529658 


-3.071799 


3.3675705 


-1.874718 


0.416538 



and triplet, singlet cases. Our fitted curves agree very 
well with the results of Jones. Somewhere, our calcu- 
lated two-term error corrected results cut and touch the 
curves of benchmark data. It is very important for HPW 
approach that our two-term error corrected results are 



free from any kind of scale, except triplet SDCS results 
at 27.2eV which have been scaled by a factor of 0.03. In 
earlier, DDCS and SDCS calculations, we had multiplied 
our data by a suitable factor for lowering the magnitude 
and compared with experimental values. 



TABLE IL CoefBcients for fitted curve at various impact en- 
ergies and singlet case. 



Energy 


a 


b 


c 


d 


27.2eV 


0.0405 


0.00567 


0.20568 


0.25395 


40.8eV 


0.018326 


0.001411 


0.049555 


0.864263 


54.4eV 


0.00373 


-0.000005 


0.00694 


0.74947 



VII. CONCLUSION 

In HPW approach, we calculate the radial part of final 
state wave function numerically, which is very crucial. 
For evaluate appropriate cross-section results, it is essen- 
tial to compute the radial part of final wave function 
very precisely. The condition of convergence depends 
on several parameters for full electron-hydrogen prob- 
lem. Model calculation is a simplification of the exact 
problem, a few parameters involved here. Currently, we 



tested the dependence of the calculation of radial wave 
function on the step length. In the figures presented in 
the paper, we have seen that with the reduction of step 
length, calculated SDCS results were better (smooth and 
less magnitude) . By using the values of T^„, elements for 
three different step lengths, we have been able to calcu- 
late two-term error corrected SDCS results. Comparison 
of two-term error corrected SDCS results with that of 
for three different step lengths shows that our endeavor 
to calculate error corrected results has been fruitful. Our 
computed error corrected results are less satisfactory, still 
there are some oscillations. Although the magnitude of 
our evaluated results quite relevant except for extreme 
asymptotic energy region. The main difficulty is that 
when we diminish the step length, the number of mash 
points is increased so there is a limitation of digital ma- 
nipulation. The fitted curves corresponding equipped 
error corrected SDCS data agree excellently with the 
benchmark results of Jones et al. [IQ] . 
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